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Abstract 

We study the electron propagator in quantum electrodynamics in lower dimensions. In the case 
of free electrons, it is well known that the propagator in momentum space takes the simple form 
Spip) = 1 / ("f ■ p — m) . In the presence of external electromagnetic fields, electron asymptotic states 
are no longer plane-waves, and hence the propagator in the basis of momentum eigenstates has a 
more intricate form. Nevertheless, in the basis of the eigenfunctions of the operator (7 • 11)^, where 
is the canonical momentum operator, it acquires the free form Sf{p) = 1/(7 -p — rn) where 
depends on the dynamical quantum numbers. We construct the electron propagator in the 
basis of the (7 • 11)^ eigenfunctions. In the (2+l)-dimensional case, we obtain it in an irreducible 
representation of the Clifford algebra incorporating to all orders the effects of a magnetic field of 
arbitrary spatial shape pointing perpendicularly to the plane of motion of the electrons. Such an 
exercise is of relevance in graphene in the massless limit. The specific examples considered include 
the uniform magnetic field and the exponentially damped static magnetic field. We further consider 
the electron propagator for the massive Schwinger model incorporating the effects of a constant 
electric field to all orders within this framework. 
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I. INTRODUCTION 



A powerful technique to solve wave equations in quantum theory is the Green's func- 
tion method. The evolution of the wave functions is encoded in the propagators, which in 
plain words are the inverses of the differential wave operators. For free particles, which are 
described by plane-waves, these propagators acquire simple diagonal forms in momentum 
space due to the isotropy of the space. However, when the isotropy is lost, particles cease 
to be described as plane-waves and their corresponding propagators have more intricate 
forms in the basis of momentum eigenstates, sometimes at the point of being intractable for 
actual calculations. An example of such systems is a gas of electrons subject to an intense 
magnetic field in a fixed direction. The non-relativistic quantum problem for this system, 
known as the Landau problem,- reveals interesting features of the dynamics of electrons in 
external magnetic fields. On one hand, the Lorentz force lacks of component along the di- 
rection of the field. This makes a clear separation between the parallel (free) and transverse 
(dynamical) components of the trajectories of the electrons with respect to the field lines. 
On the other hand, in the transverse plane, the trajectories are confined around the field 
lines. In such a case, the energy levels, called Landau levels, develop a discrete spectrum. 
These features combined translate into a rather intricate form of the electron propagator- 
as compared with its free counterpart. For relativistic electrons, the problem does not get 
any simpler. Yet, the task of unveiling these propagators for systems just as complicated is 
worthwhile. 

There have been a number of strategies developed in the past to derive the propagator 
for relativistic electrons in background electromagnetic fields, perhaps the best known being 
the Fock-Schwinger "proper time" method,-!^ which was originally developed in relativis- 
tic quantum field theory (see Ref. [5] for a brief discussion of the method), but that has 
been applied to a number of problems in ordinary non-relativistic quantum mechanics.-i^ 
However, there have been interesting alternatives to the Fock-Schwinger method to incorpo- 
rate the effects of the background fields into the propagators, like the second quantization 
of solutions to the Dirac equation in the background fields,- path integral methods,- and 
the Ritus eigenfunction method.- In this article, we derive the massive electron propaga- 
tor in background electromagnetic fields in lower space-time dimensions through the later 
method, which is based upon the diagonalization of the Dirac operator on the basis of the 
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eigenfunctions of the operator (7-11) , with 11^ being the canonical momentum operator. 

We begin this article by considering the ci-dimensional free electron propagator in Sect.HTl 
The details of the Ritus method are presented in Sect. IIIII For the sake of simplicity, 
we present the study of the electron propagator in (2+l)-dimensional QED, where a third 
spatial dimension is suppressed. This is not a mere theoretical simplification, and we explain 
ourselves: back in time, some twenty years ago,— it was shown that the low-energy effective 
theory of graphene in a tight-binding approach is the theory of two species of massless 
Dirac electrons in a (2-|-l)-dimensional Minkowski space-time, each on a different irreducible 
representation of the Clifford algebra. The isolation of graphene samples^ in 2004 and 2005, 
has given rise to the new paradigm of relativistic condensed matter,— bringing a new boost, 
both theoretical and experimental, to the matching of common interests of the condensed 
matter and high energy physics communities. Thus, the massless limit of our findings is 
of direct relevance in this subject.— We assume the electrons moving in a magnetic field 
alone pointing perpendicularly to their plane of motion. We first develop the general case 
and then, in Sect. IIVI we present a couple of examples: the motion of electrons in a uniform 
magnetic field, which is a canonical example to present the Ritus method,— and the case of 
a static magnetic field which decays exponentially along the x-axis. The effects of such a 
magnetic field to all orders can hardly be considered with other approaches, but it can be 
straightforwardly studied within the Ritus method. 

Electric field configurations can also be studied in this framework. With this goal in 
mind, in Sect.|V]we analyze the electron propagator in a constant electric field in the massive 
Schwinger model^, that is, quantum electrodynamics in (l+l)-dimensions. The Schwinger 
model exhibits spontaneous symmetry breaking of the U{1) symmetry, as well as confinement 
of electrons, thus is a favorite toy model to study some aspects of quantum chromodynamics. 
These properties can be studied from the electron propagator, which acquires a free form 
when expressed in the basis of the eigenfunctions of the operator (7-11)^. We settle the 
basis for such studies in this article. In Sect. VI we present a discussion and final remarks. 
At the end, we propose a project for the interested reader, in which the Ritus method has 
to be applied in order to obtain the electron propagator in a different situation than the 
considered here. 
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II. FREE DIRAC PROPAGATOR 

We start our discussion from the equation of motion for free electrons, namely, the Dirac 
equationi^ 



where the 7'^-matrices fulfill the Clifford algebra {'y'^,^"} = 2g'^'^, where g'^'^ = 
diag(l, — 1, . . . , — 1), depending upon the number d of space-time dimensions considered. 

In order to solve the Dirac equation, we introduce the Green's function or propagator for 
this wave equation. Propagators are important because they allow to have a visual descrip- 
tion of interaction processes out of which we desire to calculate several dynamical quantities, 
like decay rates, scattering cross sections, etcetera. In scattering problems, for instance, we 
consider a wave packet which in the remote past represents a particle approaching a poten- 
tial, and we wonder how the wave will look like in the remote future, after the interaction 
with the potential. For the sake of illustration, let us visualize this process in ordinary 
quantum mechanics. Imagine that we know the wave function ip{x, t) at one particular time 
t. At this moment, each point of space x may be viewed as a source of spherical waves 
propagating outwards in such a way that the strength of the wave amplitude arriving at 
point x' at a later time t' will be proportional to the original wave amplitude ip{x,t). This 
is Huygens' principle at work. If we denote the constant of proportionality by iG{x', t'; x, t), 
the total wave arriving at the point x' at time t' will be 



Here, G{x' ,t'; x,t) is precisely the Green's function or propagator. Thus, the knowledge of 
the Green's function enables us to construct the physical state which develops in time from 
any given initial state. In a theory without interactions the propagator is referred to as the 
free propagator, G free{x' ,t'; x,t). However if a potential V{xi,ti) is turned on at a time ti 
at the point Xi for an interval Ati, then the wave function and propagator will be modified 
as 



{i'j ■ d — m)ip{x) = 0, 



(1) 




(2) 



^{x\t') 
G{x',t';x,t) 




'Xl Gfree{x',t';Xi,ti)V{Xi,ti)^free{Xi,ti)Ati , 
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where the first term in the last expression represents the propagation of a free particle from 
{x, t) to (x', t') and the second term, read from right to left, represents free propagation from 
(x, t) to followed by a scattering at and free propagation from to 

{x',t'). The same idea might be straightforwardly translated to the relativistic case. 
In order to find the Green's function G{x,x') for the Dirac equation ([T]), we consider 

{i'y ■ d — m)G{x,x') = 5{x — x'). (4) 

Since the Green's function is translationally invariant, G{x,x') = G{x — x'), we can find its 
representation in momentum space through its Fourier transformation, Sf{p), defined as 

the asterisk denoting complex conjugation. Upon acting with the Dirac operator on G{x, x'), 
we have that 



{z^.d-m)G{x,x') = (z^.d-m) j -0-^e-'^-^SF{p){e-'^-^'y 



d'^p 



Here, in the last line we have made use of the representation of the ci— dimensional (5-function 
in terms of plane-waves. Hence, we see that 

Sf{p) = . (7) 

7 ■ p — m 

This is the free electron propagator in momentum space, suitable for calculations in quantum 
electrodynamics. 

In the next section we consider the motion of electrons in magnetic fields. A method shall 
be presented to render the propagator similar in form to the free propagator given above. 



III. PROPAGATOR IN MAGNETIC FIELDS 

The discussion in this section takes place in a (2+l)-dimensional Minkowski space-time 
with metric g^u = g^^ = diag(l, —1, —1) with a system of natural units h = c = 1 where 
c is the planar "speed of light" (Fermi velocity), which is some two orders of magnitude 
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smaller than c.—^ As we mentioned before, this is not a mere theoretical simplification. On 
the contrary, we aim to describe the low-energy theory of graphene in the massless limit. 
We start from the free Dirac equation shown in Eq. ([T]), for which we only need three 7^^- 
matrices. The lowest dimensional representation of these matrices is 2 x 2, and hence we 
can choose them to be proportional to the Pauli matrices, 

7° = (^s, 7^ = icTi, 7^ = ^cr2 , (8) 

which verify 7'^7'^ = g^^"^ — ie^^'^'^x, with 7a = Qx^"^^- Notice that there exists a second 
inequivalent representation of the Clifford algebra, specified by the 7'^-matrices 

7° = CTs, 7^ = 7^ = -^0-2 • (9) 

In graphene, the two representations describe two different electron species in each of the 
two triangular sub-lattices of the honeycomb lattice. These are conveniently merged in a 
4-component spinor with a reducible representation of the Dirac matrices given, for example, 

as 

\ -^3 y \ -a, y 

i = 1,2. For the purposes of this article, we shall only consider the representation given 
in Eq. ([8]), but the generalization is straightforward, and will be left as an exercise for the 
interested reader. 

In a background electromagnetic field, the Dirac equation takes the form 

(7 ■ n - m) ^ = 0, (11) 

where 11^ = -|- eA^ and is the electromagnetic potential defining the external field. 
Let us consider a magnetic field alone pointing perpendicularly to the plane of motion of the 
electrons, in such a fashion that, working in a Landau-like gauge, A'^ = (0, 0, W{x)), where 
W{x) is some function such that W'{x) = dxW{x) defines the profile of the field. 

As before, we are interested in solving Eq. (ITT!) through the Green's function method. In 
the presence of an external electromagnetic field, the Green's function G{x, x') satisfies 

{{-fU)-m)G{x,x')=6{x-x). (12) 

Since (7 ■ 11) does not commute with the momentum operator, neither the wave function 
nor G{x,x') can be expanded in plane-waves, and this does not allow to have a diagonal 

7 



propagator in momentum space. The scheme we choose to deal with the external fields 
was developed by Ritus.- The crucial point is to realize that the Green's function above 
should be a combination of all scalar structures obtained by contracting the 7^-matrices, the 
canonical momentum 11^ and the electromagnetic field strength tensor F^j^ = [Il^,Ily]/e = 
d^A^ — dyA^, which are compatible with Lorentz symmetry, gauge invariance and charge 
conjugation properties of quantum electrodynamics. In our conventions, we have that 

x') = G(7^n^ a^'^F^,, {F'^Ii^f) (13) 

where a^y = i[^^,-^y\/2 and 

F, = \e,yo.F^- (14) 

is the dual field strength tensor, which in (2+l)-dimensions is simply a vector. The key 
observation is that all the above structures commute with (7 ■ 11)^, that is 

[(7 ■ ^)^ (7 ■ n)] = [(7 ■ ^)^ a^^F,,] = [(7 ■ ^)^ {F^u.f] = , (15) 

relations which are straightforward to prove and represent an excellent exercise we invite 
the reader to carry out. As a consequence of the above relations, we have that 

[{^■Uy,G{x,x')]=0. (16) 

This is the main point in the Ritus method, since we know that commutating operators have 
simultaneous eigenf unctions. For our particular problem, this statement allows us to expand 
the Green's function G{x,x') in the same basis of eigenfunctions of (7 ■ 11)^. Furthermore, 
if we perform a similarity transformation on (7-11)^ in which it acquires a diagonal form in 
momentum space, then the same transformation makes the Green's function diagonal too. 
This is a general statement with operators, and hence the Ritus method can be applied to 
derive the propagator of charged scalar— and vector— particles in a similar way. 

The similarity transformation that makes (7-11)^ diagonal in the momentum space is 
such that 

E;i(7-n)X = /i, (17) 

where Ep are the transformation matrices, I is the unit matrix and can be any real 
number. Therefore, when we apply Ep functions to the propagator, it will become diagonal 
in momentum space. 
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It is important to notice that in the fermionic case, the spin operator is reahzed in terms 
of the 7'^-matrices, and thus the Ep functions inherit its matrix form. For different charged 
particles, the spin operator is reahzed in a different ways. For example, for scalar particles, 
the Ep functions are simply scalars,— whereas in the case of charged gauge bosons, the spin 
structure is embedded in a Lorentz tensor, and therefore the Ep functions also comply a 
Lorentz tensor structure.— 

Our task now is to study the structure of the Ep matrices for the case of Dirac fermions 
in (2+l)-dimensions. The similarity transformation (fTTI) can be written as 

(7 ■ U)% = p% , (18) 

which is an eigenvalue equation for the matrices Ep, which we shall refer to as the Ritus 
eigenfunctions. Explicitly, we have that 

(7 ■ n)' = 

= [7^ 7l + \ {Y, 7^) [n„ n.] + 1 {n„ n.}) 

= n2 + l^'-'-F,, . (19) 

For the configuration of the field we are considering, the only non- vanishing elements of the 
field strength tensor are F12 = — -F21 = W{x). Moreover, in the representation of the 
Dirac matrices, cr^^ = (T3, and hence the Ep functions satisfy 

(n^ + easW'ix)) Ep = Ep . (20) 

In order to solve this equation, the key observation comes from the fact that the operator 
(7 ■ n) satisfies 

[(7 ■ n),tdt] = [(7 ■ n) , -tdy] = [(7 . n) ,7^] = , (21) 

with H = — (7 ■ n)^ + IIq. Thus, these operators share a complete set of eigenfunctions, 
namely, Ep. Let the eigenvalues of these operators be 

idtEp = poEp , idyEp = -psEp , HEp = kEp . (22) 

These eigenvalues label the solutions to the massless Dirac equation in the background field. 
Furthermore, from Eq. (IT5]) . we have that = — k. Hence, the Ep functions verify 

{-Ul -Ul + easW'{x)) Ep = -kEp . (23) 
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The first two terms of tfie operator on tlie l.li.s. of this equation act on the orbital degrees 
of freedom of the eigenf unctions Ep, whereas the last term acts only in its spin degrees of 
freedom. Hence we can make the ansatz 

Ep = Ep^„ , (24) 

where uj^ is the matrix of eigenvectors of with eigenvalues a = ±1, respectively. Further- 
more, since the above operator equation is independent of t and we can look for solutions 
of the form 

■^p,a — Nf^e ^^^^ F}^ p^ „ , (25) 

with N(j being the corresponding normalization constant. Substituting the ansatze 
and fl25|) into Eq. fl23l) . we arrive at 

[dl - {-P2 + eW{x)Y + eaW'ix)) Ffc,^,,. = -kF,,p,,^ , (26) 

which is the Pauli Hamiltonian with the constrained vector potential, mass m = 1/2 and 
gyromagnetic factor g = 2, and turns out to be supersymmetric in the Quantum Mechanical 
sense (SUSY-QM).i^ 

From the solutions to the above equation, we can construct the Ritus eigenfunctions Ep 

as 

E,^f^-'=' ° V (27, 

where p = {po,P2, k) and z = {t, x, y). Being a complete set, the eigenfunctions Ep given in 
Eq. ([27]), satisfy 

j dz%>{z)'Ep{z) = 15{p-p') , 

J dp Ep{z)Ep{z') = IS{z - z') , (28) 

with Ep[z) = 7°Ep(z)7° and I is the 2 x 2 unit matrix. 

The three-momentum p^ = (po, 0, ^/k) plays an important role in the method. Its defini- 
tion involves the dynamical quantum numbers Pq and k, but not p2, which merely fixes the 
origin of the x coordinate. This means that in the Ritus method the propagator is written 
only in terms of the eigenvalues of the dynamical operators commuting with (7 ■ 11) . Notice 
that this special three-vector verifies = p^ — k = p^, and it is defined through the relation 

(7-n)Ep = Ep(7-p) . (29) 
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Indeed, we can prove this relation since, in matrix form 

(7-n)Ep = 



[-d^ + {-idy - eW)]Ep^i -idfEp^^i 

PoEp^i [~d^ - (p2 - eW)]Ep^^i 

[-d^ + ip2 - eW)]Ep^i -PoEp,^i 



(30) 



and 



E,h.p)H 'f" \ (31) 

VkEp^^i -po-Ep-i 

Here, we have made use of the properties of the Ep functions, Eq. fl22l) . We immediately 
observe that the relation holds for the diagonal components, while the off-diagonal elements 
lead to the system of equations 

l-d, - {P2 - eW)]Ep^^i = -VkEp^i , 
i-d. + iP2 - eW)]Ep,i = VkEp^^, . (32) 

If we act on the left of the second of the above expressions with the differential operator on 
the l.h.s. of the first one, we have 

l-d, - {P2 - eW)] [-9, + (p2 - eW)] Ep^, = [-d, - [p^ - eW)] VkEp,., 
[dl - {p2 - eWf + ed^W - eWd,] Ep,^ = -kEp^^ 

[dl - {p2 - eWf + eW] Ep^^ = -kEp^^ (33) 



which is nothing but Eq. ( l26l) for Ep^i. Analogously, acting on the left of Eq. (13211 with the 
differential operator on the l.h.s. of Eq. (!32|) . we obtain 

[dl + {p2 - eWf - eW']Ep^^i = -kEp^^,, (34) 

which is again Eq. ( l26l) . now for Thus the relation (!29l) is also valid for the off-diagonal 

elements. 

With the Ep functions, we can consider the Green's function method to obtain the prop- 
agator in the presence of the field. We start from Eq. f|T2|) and re-label {x,x') {z,z'). 
Then we define the Green's function in momentum space as 

G{z,z')= J dpEp{z)SF{pWp{z') . (35) 
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We would like to stress that the integration might as well represent a sum, depending 
upon the continuous or discrete nature of the components of the momentum. Applying the 
Dirac operator (7 ■ 11 — m) to G{z, z'), we have that 

{-fU-m)G{z,z') = (7-n-m) j dp Ep{z)Sf{p%{z') 

dp (7 ■ n - m)Ep{z)SF{p)Ep{z') 
dp Ep(2;)(7 ■ p - m)SFip)Ep{z') 
= j dp ¥.p{z)%{z') , (36) 

where in the last line we have used the representation of the 5-function in the Ep basis. 
Hence we notice that the Ritus propagator takes the form of a free propagator, namely, 

Sf{p) = , (37) 

7 ■ p — m 

in this basis, with p^ defined through Eq. ( l29l) . 

On physical grounds, the Ep functions correspond to the asymptotic states of electrons 
with momentum p in the background of the external field under consideration. With the help 
of these functions and the property fl29l) . we can find the solutions of the Dirac equation ffTTj) 
in a straightforward manner. To this end, we propose \E' = EpWp, then 

(7 ■ n — m)'EpUp = Ep(7 ■ p — m)up = 0, (38) 

and thus we see that Up is simply a free spinor describing an electron with momentum 
p. Notice that with this form of \E', the information concerning the interaction with the 
background magnetic field has been factorized into the Ep functions and throughout the 
p dependence of Up. Several relevant physical observables can then be found immediately, 
such as the probability density, the transmission and reflection coefficients between magnetic 
domains, and the density matrix, which are all useful, for example, in graphene applications 
as those we discuss below, where the Ritus method plays useful. 
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FIG. 1: The sketch of the field hnes of a uniform magnetic field. 



IV. EXAMPLES 

A. Uniform Magnetic Field 

To have an explicit taste of the power of Ritus method, let us first consider the case 
of a uniform magnetic field.— This corresponds to the choice W{x) = Bqx and the field 
lines are sketched in Fig. [TJ To simplify the calculations, we rename the quantum number 
k — >■ 2\eBQ\k in Eq. fl22|) . In this case, Eq. fl26l) simplifies to 



[dl - {-P2 + eBoxf + aeBo] Fk^p.^^) = -2\eBo\kFk^p,^^{x) . (39) 
Let 7] = ^/2\eBo\{x — Then the above expression acquires the form 



+ k + -sgn(eSo) - 



Fk,p,Av) = , (40) 



_9r/2 ■ 2 ° ' 4 

that is, the equation for a quantum harmonic oscillator, with center of oscillation in Xq = 
P2/(e-Bo) and cyclotron frequency Wc = 2eBo. Thus, the normalized functions Ep^^ acquire 
the form 

where 

Dr,{x) = 2-"/2e-V4^^ (^x/Vtj (42) 
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FIG. 2: The Ep functions, Eq. (j41|) . in arbitrary units along the dynamical direction for various 
Landau levels. The scale of the graphic is set by cBq = 1 and p2 = 0. 

is the paraboUc cyUnder function of order 



and Hn{x) are the Hermite's polynomials. Expectedly, the uniform magnetic field renders 
the {n — l)-th state with spin down with the same energy of the n-th state with spin up. 
In Fig. [2], we depicted the functions Ep = Ep^i along the dynamical direction for various 
Landau levels. 

As we mentioned before, the Ritus eigenf unctions correspond to the asymptotic states of 
electrons in the presence of a magnetic field. Therefore, physical observables like probability 
densities are linear combinations of l-Epp. These functions have the profile shown Fig. [31 
The solid curve enveloping these solutions corresponds to the potential 



where W = eW — p2 is referred to as the superpotential in the SUSY-QM literature.— 
B. Exponential Magnetic Field 

There are many problems relating electrons in non-uniform magnetic fields of relevance 
in graphene. In particular, it has been established the possibility to confine quasiparticles in 
magnetic barriers.— This could be feasible creating spatially inhomogeneous, but constant in 
time, magnetic fields depositing ferromagnetic layers over the substrate of a graphene sample 
layer.— The physical properties of graphene make it a promising novel material to control 




(43) 




(44) 
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FIG. 3: The l-EpP function from Eq. (j4ip in arbitrary units along the dynamical direction for 
various Landau levels. The curve y = — 1 corresponds to the potential (j44p for this field 
configuration. The scale of the graphic is set by cBq = 1 and p2 = 0. 




FIG. 4: The sketch of the field lines of an exponentially decaying static magnetic field. 

the transport properties in nanodevices. It has been considered to be used in electronics and 
spintronics applications, like in single-electron transistors,— in the so called magnetic edge 
states,— which may play an important role in the spin-polarized currents along magnetic 
domains, and in quantum dots and antidots magnetically confined. Moreover, the quantum 
Hall effect in graphene has been observed at room temperature,— evidence which confirms 
the great potential of graphene as the material to be used in carbon-based electronic devices. 

Here we study the electron propagator in a background static magnetic field which has 
an exponentially decaying spatial profile along one direction, described through the function 
W{x) = — -Bo exp{— ax}. The field lines of such magnetic configuration are sketched in 
Fig. m In this case, Eq. fl26|) simplifies to 

[dl - {-p2 - e5oexp{-Q!x})^ + aaeBo exp{-ax}] Ffc_p2,^(x) = -kFk,p^^^{x) . (45) 
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FIG. 5: The Ep functions from Eq. (j47p in arbitrary units along the dynamical direction for various 
values of s at fixed n = 0. The scale of the graphic is set by cBq = a = 1. 

Let Q = (2ei?o/tt) exp{—ax} and s = —p2/(y, then, the above expression is equivalent to 
This equation has the normalized solutions Ep^^j given as 

P _ 1 ( 2an!(5-n) V''^ jp^t+jp^y /2 U-n)j2{s-n)( 

"'^ ~ 271 \T{2s-n + l)) e Q {q), 

where L^^i^) are the associate Laguerre polynomials with 

n = .-yiA7^. (48) 

The quantum number n is the principal quantum number, whereas s a center of oscillation 
weighted by the damping factor a. The solutions Ep from Eq. P7|) are sketched in Figs. [MZl 
for n = 0, 1 and 2 for various values of s, whereas in Fig. [8] we show l-EpP for various values 
of n at fixed s = 8. Notice that in this case the potential (jH]) also envelops the squares of 
the solutions. 

This concludes our presentation of the Ritus method in magnetic fields. In the next 
section we consider the Schwinger model for electric fields within this approach. 

V. PROPAGATOR IN ELECTRIC FIELDS 

The Dirac equation in the massive Schwinger model^ takes the form shown in Eq. (fTTIl . 
but only two 7'^-matrices are involved. In this section, we work in a (l+l)-dimensional 
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Fk,p,A0)=O- (46) 



Epix), n=l 




X 



FIG. 6: The Ep functions from Eq. (j47p in arbitrary units along the dynamical direction for various 
values of s at fixed n = 1. The scale of the graphic is set by cBq = a = 1. 



Epix), n=2 




X 



FIG. 7: The Ep functions from Eq. (|47p in arbitrary units along the dynamical direction for various 
values of s at fixed n = 2. The scale of the graphic is set by cBq = a = 1. 

Minkowski space with metric g'^'^ = g^y = diag(l, —1), again in natural units. We choose 
the 7'^-matrices as 



7o 



1 

1 



7i 



-1 

1 



(49) 



Notice that 



75 = 7o7i 



1 
-1 



o"oi • 



(50) 



We introduce a constant uniform electric field in the gauge 



A^ = {-Eox,0) 



(51) 
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FIG. 8: The \Ep\'^ functions from Eq. (07]) in arbitrary units along the dynamical direction for 



various values of n at fixed s = 8. The curve y 



~x\2 



+ e ^ corresponds to the potential 



for this field configuration. The scale of the graphic is set by cBq = a = 1. 

and then, since the only non-vanishing components of the field strength tensor are Fqi 
—FiQ = Eq, we have that 



(7 ■ UfEp = [U^ - lej.Eo] = p^E, 
We look for solutions of the form 

Ep = Ep^^{t,x)uJa 



(52) 



e-'P'' Ep,^{x)uj^ , 



(53) 



where is the matrix of eigenvectors of 75 with eigenvalues a = ±1, respectively. Explicitly, 
the Ritus eigenf unctions for this problem are of the form 



Ep = e 



-ipot 




(54) 



where the Ep^cr functions verify 

[-{dt + leEoxY + dl- leaEo] Ep^„ = p^Ep^, . (55) 
Let us make the replacement p^ — » 2ieEQp^ as in the case of a uniform magnetic field. Then 



making the change of variables Eq iEq and p = y2e8Q {x + ip^/ (e£^o)), we get 

EpAp) = , (56) 



92 2 ^ P^' 
+2- 4 
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which is again a displaced harmonic oscillator with center of oscillation in xq = ipo/ieSo) 
and frequency Uc = 2eSo- Thus the solutions are again parabolic cylinder functions 

^P.i(P) = (27r)3/2(n + l)!i/2^"+^^^) ' 
EpMp) = ^^^.^nip) , (57) 

of order 

^ = P' + |-^- (58) 

Contrary to the case of the uniform magnetic field, the arguments of the parabolic cylinder 
functions as well as the order, are complex numbers. Nevertheless, they serve as well to 
diagonalize the propagator as 



S(x) = [ dp Ep(x) — =i EJx) , (59) 

J •y ■ p — m 

which is found through the relation (!29|) with p^ = (p, 0), where p is defined in Eq. ( 158|) . 

This illustrates the usefulness of the Ritus method. Below we present a discussion and 
the conclusions. 



VI. DISCUSSION AND CONCLUSIONS 

In this paper we have reviewed the massive electron propagator in lower than (3+1) 
space-time dimensions. We have obtained the free electron propagator through the Green's 
function method in arbitrary dimensions and has, as we know, a diagonal form in momentum 
space, namely Sf{p) = 1/(7 ■ p — m). This is obviously the case because the free Dirac 
operator commutes with the momentum operator and thus one can take simply a Fourier 
transform, that is to say, we can use the eigenfunctions of the momentum operator to 
diagonalize the Green's function. 

We have then considered the electron propagator in (2+l)-dimensions in an external 
magnetic field pointing perpendicularly to the plane of motion of the electron and whose 
spatial shape we considered described by some arbitrary function W{x). This, to the best 
of our knowledge, has been presented for the first time in literature. We selected to work 
within the Ritus method^ framework, and construct explicitly the eigenfunctions Ep of the 



19 



(7 ■ n)^ operator in the general case. The basic idea behind this method is the following: 
since the Dirac operator (7 ■ 11) does not commute with the momentum operator, we cannot 
simply take the Fourier transform of the Green's function and have a diagonal propagator in 
momentum space. However since (7 ■ 11) commutes with idt-, —idy and, more importantly, 
with (7 ■ n)^, we can use their common eigenfunctions, Ep, to diagonalize the propagator, 
which turns out to have the form of a free propagator with momentum depending upon the 
dynamical quantum numbers. We have specialized our findings to the case of a uniform 
magnetic field where the Ep functions are described in terms of parabolic cylinder functions. 
As a second example, we considered an exponentially damped static magnetic field. Here 
the Ritus eigenfunctions are written in terms of associate Legendre functions. On both these 
cases, the massless versions of the propagators are of direct relevance in graphene. 

We also considered the propagator in (l+l)-dimensions in a uniform electric field. Our 
findings in this case are similar to the magnetic field case in the plane, except that the 
eigenfunctions have complex arguments and order. The important lesson here is that the 
method works in arbitrary dimensions and different field configurations. In each case, the 
propagator continues to have a diagonal form, although the Ep functions change from one 
configuration to the other. 

In conclusion, the Ritus method offers an alternative to study the electron propagator 
in the presence of external fields. Contrary to the "proper time" method^i^, the propagator 
here acquires its free form, and the effects of the field are reflected in the eigenfunctions 
of the operator (7 ■ H)^, namely Ep. The method works as well in the case of scalar— and 
vector charged particles.— The idea behind it is familiar in quantum mechanics, where the 
diagonalization of a given operator is carried out with the aid of eigenfunctions of a second 
commuting operator. That makes Ritus method suitable to study the propagator for more 
complex configurations of external fields. 

Project 

For the interested reader, we propose the following project: Work out the derivation of 
the propagator in (2+l)-dimensions in a magnetic field making use of the second irreducible 
representation of the Dirac matrices, Eq.Q, and show that in this case, the Ep functions 
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have the form 



E,= |^"-''=> ° |, (60) 
E.i(2) 
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